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Abstract: Models of natural supersymmetry seek to solve the little hierarchy problem by 
positing a spectrum of light higgsinos < 200 — 300 GeV and light top squarks < 600 GeV 
along with very heavy squarks and TeV-scale gluinos. Such models have low electroweak 
fine-tuning and satisfy the LHC constraints. However, in the context of the MSSM, they 
predict too low a value of m^, are frequently in conflict with the measured 6—7-57 branching 
fraction and the relic density of thermally produced higgsino-like WIMPs falls well below 
dark matter (DM) measurements. We propose a framework dubbed radiative natural SUSY 
(RNS) which can be realized within the MSSM (avoiding the addition of extra exotic 
matter) and which maintains features such as gauge coupling unification and radiative 
electroweak symmetry breaking. The RNS model can be generated from SUSY GUT 
type models with non- universal Higgs masses (NUHM). Allowing for high scale soft SUSY 
breaking Higgs mass mn u > m o leads to automatic cancellations during renormalization 
group (RG) running, and to radiatively-induced low fine-tuning at the electroweak scale. 
Coupled with large mixing in the top squark sector, RNS allows for fine-tuning at the 3-10% 
level with TeV-scale top squarks and a 125 GeV light Higgs scalar h. The model allows for 
at least a partial solution to the SUSY flavor, CP and gravitino problems since first/second 
generation scalars (and the gravitino) may exist in the 10-30 TeV regime. We outline some 
possible signatures for RNS at the LHC such as the appearance of low invariant mass 
opposite-sign isolated dileptons from gluino cascade decays. The smoking gun signature 
for RNS is the appearance of light higgsinos at a linear e + e~ collider. If the strong CP 
problem is solved by the Peccei-Quinn mechanism, then RNS naturally accommodates 
mixed axion-higgsino cold dark matter, where the light higgsino-like WIMPS - which in 
this case make up only a fraction of the measured relic abundance - should be detectable 
at upcoming WIMP detectors. 



1. Introduction 

The recent discovery by Atlas and CMS of a Higgs-like resonance at the CERN LHC[1, 2] 
adds credence to supersymmetric models (SUSY) of particle physics in that the mass 
value nth — 125 GeV falls squarely within the narrow window predicted by the Minimal 
Supersymmetric Standard Model (MSSM): my L ~ 115 — 135 GeV[3]. At the same time, the 
lack of a SUSY signal at LHC7 and LHC8 implies m g > 1.4 TeV (for m g ~ rrig) and m g > 
0.9 TeV (for m g <C m g )[A, 5]. While weak scale SUSY[6] provides a solution[7] to the gauge 
hierarchy problem via the cancellation of quadratic divergences, the apparently multi- 
TeV sparticle masses required by LHC searches seemingly exacerbates the little hierarchy 
problem (LHP): 

• how do multi-TeV values of SUSY model parameters conspire to yield a Z-boson 
(Higgs boson) mass of just 91.2 (125) GeV? 

Models of natural supersymmetry[8] address the LHP by positing a spectrum of light 
higgsinos < 200 GeV and light top- and bottom-squarks with g < 600 GeV along 
with very heavy first/second generation squarks and TeV-scale gluinos[9, 10, 11]. Such a 
spectrum allows for low electroweak fine-tuning (EWFT) while at the same time keeping 
sparticles safely beyond LHC search limits. Because third generation scalars are in the few 
hundred GeV range, the radiative corrections to rrih, which increase logarithmically with 
m~ , are never very large and these models have great difficulty in accommodating a light 
SUSY Higgs scalar with mass rrih ~ 125 GeV [11, 12]. Thus, we are faced with a new 
conundrum: 

• how does one reconcile low EWFT with such a large value of rrih [13]? 
A second problem occurs in that 

• the predicted branching fraction for b — > sj decay is frequently at odds with the 
measured value due to the very light third generation squarks [11]. 

A third issue appears in that 

• the light higgsino-like WIMP particles predicted by models of natural SUSY lead to 
a thermally-generated relic density which is typically a factor 10-15 below[14, 11] the 
WMAP measured value of QcDAih 2 ~ 0.11. 

One solution to the fine-tuning/Higgs problem is to add extra fields to the theory, thus 
moving beyond the MSSM[13]. For example, adding an extra singlet as in the NMSSM 
permits a new quartic coupling in the Higgs potential thus allowing for an increased value 
of rrih [15]. Alternatively, one may add extra vector-like matter to increase rrih while main- 
taining light top squarks[16]. In the former case of the NMSSM, adding extra gauge singlets 
may lead to re-introduction of destabilizing divergences into the theory[17]. In the latter 
case, one might wonder about the ad-hoc introduction of extra weak scale matter multiplets 
and how they might have avoided detection. A third possibility, which is presented below, 
is to re-examine EWFT and to ascertain if there do indeed exist sparticle spectra within 
the MSSM that lead to rrih ~ 125 GeV while maintaining modest levels of electroweak 
fine-tuning. 
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1.1 Electroweak fine-tuning 

One way to evaluate EWFT in SUSY models is to examine the minimization condition 
from the Higgs sector scalar potential which determines the Z boson mass. (Alternatively, 
one may examine the mass formula for and arrive at similar conclusions.) Minimization 
of the one-loop effective potential Vt Tee + AV, leads to 

M 2 Z _ m 2 Hd + ^ d d - (ml u + Eg tan 2 g % 
2 " tan 2 /?-l M ' [ ' 

where E^ and E^ are radiative corrections that arise from the derivatives of AV evaluated 
at the minimum. Eq. (1.1) reduces to the familiar tree-level expression[6] for M| when 
radiative correction terms are ignored. As we will discuss in detail below, Ti^ and E^ include 
contributions, listed in the Appendix, from various particles and sparticles with sizeable 
Yukawa and/or gauge couplings to the Higgs sector. To obtain a natural value of Mz on 
the left-hand-side, one would like each term (with i = H^, H u , \x as well as E"(fc), 
E^(/c), where k denotes the various contributions to the Es that we just mentioned) on the 
right-hand-side to have an absolute value of order M|/2. Noting that all entries in (1.1) 
are defined at the weak scale, we are led to define the electroweak fine-tuning parameter 1 
by [20] 

A EW = maxi (d) /(Aff/2) , (1.2) 

where Ch u = \ — m2 u u tan 2 /3/(tan 2 f3 — 1)|/, Cn d = |m# /(tan 2 f3 — 1)|/ and C M = | — fi 2 \, 
along with analogous definitions for Csjj(fc) and C S d^. Low Aew means less fine-tuning. 
Since Cu d and C s <j^^ terms are suppressed by tan 2 f3 — 1, for even moderate tan/3 values 
this expression reduces approximately to 

^^-im 2 Hu + K)-^. (1-3) 

We see that to get low Aew we require | — rn 2 H J\ ~ M§/2 and /i 2 ~ M|/2. The question 
then arises: what is the model and can we find a set of model parameters such that 
Aew ~ 1 — 30, corresponding to better than A^ w = 3% EWFT? Note that Aew depends 
only on the weak scale parameters of the theory and hence is essentially fixed by the particle 
spectrum, independent of how superpartner masses arise. 

To understand how the underlying framework for superpartner masses may be relevant, 
consider a model with input parameters defined at some high scale A 3> AfsusY, where 
Msusy is the SUSY breaking scale ~ 1 TeV and A may be as high as Mqut or even the 
reduced Planck mass Mp. Then 

m| u (M SUSY ) = m 2 Hu (A) + Sm 2 Hu (1.4) 

1 Barbieri and Giudice[18] (and, even earlier, Ellis et aZ.[19]) define a fine tuning measure Abg = 
max\(a,i/Mz)dMz/dai\ for input parameters a%. Our definition coincides with theirs when M§ depends 
linearly on input parameters (such as /i 2 , m 2 Hu or m 2 Hd using electroweak scale parameters) but differs when 
the parameter dependence is non-linear. For electroweak scale parameters, the non-linear dependence only 
occurs in the radiative correction terms E£J and and in tan/3. 
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where 

*»k - S <"4 + < + A i) "* (^M • (1-5) 



Requiring 8m 2 H < Ax ^ then leads for m^ = 125 GeV to 

sin/3 Aog^V V VA^ 1/2 



rn|+m|< 600 GeV-^f^v j (±) ^ d-d 



where Rt = At/ ^Jm~ + m~ . Taking A = 10 and A as low as 20 TeV corresponds to[8, 9, 10] 
• H < 200 GeV, 



m lo m k <600 GeV, 



mj < 1.5 - 2 TeV. 



The last of these conditions arises because the squark radiative corrections 5m~ ~ (2gj /3n )m~ x 
log A. Setting the log to unity and requiring <5m 2 < m? then implies m„ < 3mr , or 
m^ < 1.5 — 2 GeV for A < 10. Taking A as high as Mqjjt leads to even tighter constraints: 
m h 2 ' m b 1 ~ ^00 GeV and m^ < 600 GeV, almost certainly in violation of LHC sparticle 
search constraints. Since (degenerate) first/second generation squarks and sleptons enter 
into (1.1) only at the two loop level, these can be much heavier: beyond LHC reach and 
also possibly heavy enough to provide a (partial) decoupling solution to the SUSY flavor 
and CP problems. In gravity mediation where rriq ~ m 3 / 2 , then one also solves the cos- 
mological gravitino problem [21, 22] and in GUTs one also suppresses proton decay. Then 
we may also have 

• m- j ~ 10 - 50 TeV. 

The generic natural SUSY (NS) solution reconciles lack of a SUSY signal at LHC with 
allowing for electroweak naturalness. It also predicts that the ti t 2 and b\ may soon be 
accessible to LHC searches. New limits from direct top- and bottom-squark pair production 
searches, interpreted within the context of simplified models, have begun to bite into the 
NS parameter space[23]. Of course, if tti^, rn^ ~ m z ■> then the visible decay products 
from stop and sbottom production will be soft and difficult to see at the LHC. 

A more worrisome problem comes from the newly discovered value of the Higgs mass 
nih ~ 125 GeV. In the MSSM, one obtains[3] (assuming that the t-squarks are not very 
split), 



,2 _ nj2 „„„2 na , ^ m t 
8tt 2 m 



m\ ~ Af§ cos 2 2/3 + 



In -4 + -L 1 



771? J7l~ \ 12m? 

t \ t . 



(1.7) 



where Xt = At — \i cot ft and m~ ~ mQ 3 mjj 3 . For a given m~, this expression is maximal for 
large mixing in the top-squark sector with X™ ax = \^6mf. With top-squark masses below 
about 500 GeV, the radiative corrections to m^ are not large enough to yield m^ ~ 125 GeV 
even with maximal mixing[13, 11]. This situation has been used to argue that additional 
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multiplets beyond those of the MSSM must be present in order to raise up while 
maintaining very light third generation squarks[13]. Added to these are the two issues 
mentioned earlier: 1. the very light third generation squarks[ll] endemic to NS lead to 
a predicted branching fraction for b — > decay which is frequently much lower than the 
measured value[ll], and 2. that the relic abundance of higgsino-like WIMPs inherent in 
NS, calculated in the standard MSSM-only cosmology, is typically a factor 10-15 below 
measured values [11]. These issues have led to increasing skepticism of weak scale SUSY as 
realized in the natural SUSY incarnation described above. 

A possible resolution to the above issues associated with a NS spectrum is to simply 
invoke a SUSY particle spectrum at the weak scale (or some other nearby scale [24]), as 
in the pMSSM model [25] so that large logarithms associated with a high value of A are 
absent. In this case, A ~ MgusY and Sm^ is not enhanced by large logarithms and we 
may select parameters m 2 H ~ /s 2 ~ M§ ~ m^. Of course, heavy top squarks are needed 
to obtain the observed value of m^. While a logical possibility, this solution loses several 
attractive features of models which are valid up to scales as high as A ~ Mgut> such as 
gauge coupling unification and radiative electroweak symmetry breaking driven by a large 
top quark mass. 

Another alternative is to use Aew defined above as a fine-tuning measure even for 
models defined at the high scale. This use of weak scale parameters to define the fine- 
tuning criterion is a weaker condition since it allows for possible cancellations in (1.4). 
Indeed this is precisely what happens in what is known as the hyperbolic branch or focus 
point region (HB/FP) of mSUGRA[26]: m 2 Hu {k) + 5m 2 Hu ~ ro^ (M S usy) ~ M 2 ~ M 2 . 
The HB/FP region of mSUGRA occurs, however, only for small values of Aq/ttiq [27] and 
yields < 120 GeV, well below the Atlas/CMS measured value of ~ 125 GeV. Scans 
over parameter space show that the HB/FP region is nearly excluded if one requires both 
low \fi\ and m h ~ 123 - 127 GeV [28, 27]. 

To obtain a viable high scale model we see that we clearly need to go beyond mSUGRA. 
The small value of —m 2 H ^(MsuSY) that we require for low EWFT can be obtained in sev- 
eral ways. For instance, we could introduce non-universality of gaugino masses and adopt 
very high GUT scale values of the SU(2) gaugino mass parameters [29], or a low value of 
the SU(3) gaugino mass parameter [30]. Both choices would lead to a larger chargino to 
gluino mass ratio than in models with universal gaugino masses [31], and since charginos 
couple directly to the Higgs sector, potentially significant contributions to the radiative 
corrections for gluinos that satisfy the LHC bound. The other way of obtaining small val- 
ues of —m 2 Hu {MsusY) without undue cancellations in (1.1) is to introduce non- universality 
in the scalar Higgs sector. To facilitate our analysis, we use the two parameter non- 
universal Higgs mass (NUHM2) extension[32] of the mSUGRA model where (Mqut) 
and m 2 H {Mq^), or equivalently, the weak scale parameters \i and m^, are chosen inde- 
pendently of matter scalar mass parameters, and the model is completely specified by the 
parameter set, 

m , mi/2, Mi m A, Aq, tan/3 (NUHM2). (1.8) 
Modest electroweak fine-tuning is then obtained due to large cancellations between m 2 H (A = 
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-^gut) and Smjj . 

Along with —m 2 H ^{MsusY) ~ ^§/2, low Aew also requires ^ 2 ~ M%/2. In gravity- 
mediated SUSY breaking models- where the /i problem is solved by the Guidice-Masiero 
mechanism [33]- one expects that fi ~ A?ti 3 / 2 where A is a hidden- visible sector coupling. 
For small values of A, then we expect <C m 3 / 2 which is then significantly smaller than 
the typical soft-SUSY breaking masses. 

Later in our analysis, we will also allow for the possibility of split matter generations 
where the third generation mass parameter mo (3) is independent from the corresponding 
parameter m-o(l,2) for the first two generations. We refer to this case as the NUHM3 
model. The NUHM3 model allows for an improved decoupling solution to the SUSY flavour 
problem. Finally, we mention that a small magnitude of n is also possible in the NUHM1 
model where we take the two GUT scale Higgs mass parameters to have a common value 

which can be raised above mo until jtt^^Mstjsy) becomes comparable to M| [34]. In 
this case, m^ is of course determined. We consider this case briefly. In all cases, intra- 
generation splitting is avoided since as noted in Ref. [27], it can lead to large fine-tuning if 
scalars are very heavy. We emphasize that Aew is determined by physical sparticle masses 
and couplings so that our results can be applied to any model that yields a similar spectrum, 
irrespective of how sparticles acquire their SUSY breaking mass. 

At this point the reader may legitimately wonder about the validity of using Aew as 
a measure of fine-tuning in the NUHM2 model which is, by construction, assumed to be a 
description of physics up to energy scales as large as I/qut- The introduction of the Ahs in 
Ref. [27] was precisely to include the impact of these large logarithms - which also appear in 
(1.6) - on the fine-tuning. The use of Aew as a fine-tuning measure allows for the possibility 
of large cancellation between m 2 H d (A) and the term 8m 2 H d (which may include large 
logarithms). For instance, special regions of parameter space of some models {e.g. focus- 
point SUSY, the mixed-modulus-anomaly-mediation model for special values of the ratio 
a [35] or particular regions of parameter space of non- universal Higgs mass models) enjoy 
nearly complete cancellations between the terms with large logarithms and m 2 H d (A). In a 
more encompassing framework that includes the origin of soft SUSY breaking parameters, 
such cancellations might not only be allowed, but might be automatic [35, 36]. We note, 
however, that Ahsj as we have defined it, does not take such a cancellation into account; 
under these circumstances Aew is the appropriate fine-tuning measure to use. 

The fine-tuning measure Aew introduced in Ref's. [20] and [27] has several attractive 
features that merit consideration. 

• Model independent (within the context of models which reduce to the MSSM at the 
weak scale): Aew is essentially determined by the sparticle spectrum[27], and - 
unlike Ahs and other measures of fine-tuning - does not depend on the mechanism 
by which sparticles acquire masses. Since Aew is determined only from weak scale 
Lagrangian parameters, the phenomenological consequences which may be derived 
by requiring low Aew win apply not only for the NUHM2 model considered here, 
but also for other possibly more complete (or less complete, such as pMSSM) models 
which lead to look-alike spectra at the weak scale. 
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• Conservative: Aew captures the minimal fine-tuning that is necessary for any given 
sparticle spectrum, and so leads to the most conservative conclusions regarding fine- 
tuning considerations. 

• Measureable: Aew is m principle measurable in that it can be evaluated if the un- 
derlying weak scale parameters can be extracted from data. 

• Unambiguous: Fine-tuning measures which depend on high scale parameter choices, 
such as the Barbieri-Guidice measure Abg discussed previously, are highly sensitive 
to exactly which set of model input parameters one adopts: for example, it is well- 
known that significantly different values of Abg result depending on whether the high 
scale top- Yukawa coupling is or is not included as an input parameter [37]. There is 
no such ambiguity in the fine-tuning sensitivity as measured by both Aew an d Ahs- 

• Predictive: While Aew is less restrictive than Ahs-, it still remains highly restrictive. 
The requirement of low Aew highly disfavors models such as mSUGRA/CMSSM[27], 
while allowing for very distinct predictions from more general models such as NUHM2. 

• Falsifiable: The most important prediction from requiring low Aew is that \^\ cannot 
be too far removed from Mz- This implies the existence of light higgsinos ~ 100 — 300 
GeV which are hard to see at hadron colliders, but which are easily detected at a 
linear e + e~ collider with yfs > 2|/i|. If no higgsinos appear at ILC1000, then the idea 
of electroweak naturalness in SUSY models is dead. 

• Simple to calculate: Aew is extremely simple to encode in sparticle mass spectrum 
programs, even if one adopts models with very large numbers of input parameters. 

To illustrate how a low value of m H (Mstjsy) is obtained, in Fig. 1 we show the running 
of various SUSY parameters versus the renormalization scale Q for the RNS2 benchmark 
point from Ref. [20]. The RNS2 point has parameters iuq = 7025 GeV, m^j = 568.3 GeV, 
A = -11426.6 GeV, tan/3 = 8.55 with /x = 150 GeV and m A = 1000 GeV. The gaugino 
and matter scalar mass parameters evolve from m^ 2 and too to their weak scale values, 
resulting in a pattern of masses very similar to that in mSUGRA. The parameter [i hardly 
evolves, and for such a low value of tan /3, ni Hd also suffers little evolution. Of most 
interest to us here is the RG evolution of m H . As is well known, the SUSY breaking 
parameters to,q 3 , ?n^ 3 and rn? H ^ of the scalar fields that couple via the large top quark 
Yukawa coupling are driven down with reducing values of the scale Q. The reduction is 
the greatest for m H which, in fact, is driven negative, triggering the radiative breakdown 
of electroweak symmetry [38]. We see from the figure that the weak scale value of —m H 
has a magnitude ~ M| , and is much smaller than the weak scale value of other mass 
parameters. This is not an accident because the NUHM2 model provides us the flexibility 
to adjust the GUT scale value of m Hu so that it barely runs to negative values at the weak 
scale. Since m H is driven radiatively to ~ — M| at the weak scale, this scenario has been 
dubbed Radiative Natural SUSY, or RNS for short. 

In Fig. 2 we scan over parameter space of the NUHM2 model - while enforcing = 
125 ± 2 GeV and LHC sparticle mass limits - and plot the value of Ahs versus the high 
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Figure 1: Evolution of SSB parameters from Mqut to M wea k for the RNS2 benchmark point 
taken from in Ref. [20] whose parameters are given in the text. The graph extends to values below 
Q 2 = m^m^ where the Higgs mass parameters are extracted. 

scale matter scalar mass parameter uiq. We see that the smallest value of Ahs is 10 3 for 
the lowest values of mo allowed. This is because of the large logarithms that we discussed 
above. The NUHM2 would be fine-tuned to at least 0.1%, and usually even higher fine- 
tuning is necessary. We will see below that much smaller values of Aew are possible in 
some parts of parameter space. An underlying high-scale theory that automatically leads 
to an NUHM2-like spectrum in this parameter region would then not be so fine-tuned. In 
contrast, there are no analogous regions of mSUGRA/CMSSM parameters for which we 
have shown that Aew ^ 100 [27]. 

In the remainder of this paper, we will explore what parameter choices lead to low val- 
ues of A E w- While A E w seems bounded from below by about 100 in mSUGRA/CMSSM[27], 
we will find that Aew as l° w as ~ 10 can be obtained in NUHM2. In addition, requiring 
low Aew 30 places strong restrictions on the allowed sparticle mass spectra, leading to 
distinctive predictions for collider and dark matter searches. 

1.2 Radiative natural supersymmetry 

Motivated by the possibility of cancellations occurring in mjj^Msusy), we return to the 
EWSB minimization condition (1.1) which was introduced earlier, and examine more care- 
fully the radiative corrections embodied in and that we have not discussed up to 
now. These affect the minimization condition in an important way when m^(msuSY) arid 
/j? are much smaller than the scale of other weak scale SUSY breaking parameters. At the 
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Figure 2: Plot of Ahs versus mo from a scan over NUHM2 model parameters, accepting only 
points where = 125 ± 2 GeV and which obey LHC sparticle mass constraints. 



one-loop level, £^ contains the contributions [39, 40] ££(*i, 2 ), £"(61,2), ^(fi >2 ), ££(Wi >2 ), 
E«(2i_ 4 ), ££(/i, H), ££(#*), EJKW*), and ££(*). contains similar terms along 

with S^(6) and S^(r) while E^(t) = 0. The complete set of one-loop contributions to these 
is listed in the Appendix. There are additional contributions from first/second generation 
sparticles from their D-term couplings to Higgs scalars. If these squarks, and separately 
sleptons, are degenerate then these contributions cancel within each generation because the 
sum of weak isospins/hypercharges of squarks/sleptons total to zero[27]. In the parameter 



space region where RNS is realized, 



i.e. 



where 



(Msusy) ~ /^ 2 



M|, the radiative 



correction terms from S" may give the largest contributions to Aew- 

The largest of the E" terms almost always come from top squarks for which we find, 



S u(*l,2 

= (m 



16vr 2 



F(m~ ) 
v 11,2' 



ft 



2 ^ tf A t 

9z^ 



83K2 " h w )A 



2 



2 

ra- 
ti 



(1.9) 



where A t = {mf L - m| )/2 + Mj cos2/3(| - %x w ), 9 Z = (sf + 9 )/8, ^VK = sin^ W 
and F(m 2 ) = m? (log(m 2 /Q 2 ) — l), with Q 2 = m^m^ . In Ref. [20], it is shown that for 
the case of the t\ contribution, as \A t \ gets large there is a suppression of £„(ii) due to 
a cancellation between terms in the square brackets of Eq. (1.9). The i 2 contribution is 
suppressed if there is a sizeable splitting between and due to a large cancellation 
within F(m~ ) because log(m| / 'Q 2 ) = log (m^/m^) ~ 1. The large |^4t| values suppress 
both top squark contributions to E^, and at the same time lift up the value of mh, which 
is near maximal for large negative A%. Combining all effects, one sees that the same 
mechanism responsible for boosting the value of rrih into accord with LHC measurements 
can also suppress the E" contributions to EWFT, leading to a model with low EWFT. 

To display the quality of EWFT explicitly, we show in Fig. 3a the various signed 
contributions to M§/2 that enter Eq. (1.1) for the RNS2 point from Fig. 1 and Ref. [20]. 
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In this figure, we label these signed contributions by Q where i = H u , Hd, (J,, EJJ, S^. The 
largest contributions come from Csjj ~ 0.04 TeV 2 and Cjj u ~ —0.03 TeV 2 . In frame b), we 
show these same quantities for the mSUGRA model (where fi and tua are outputs instead of 
input parameters). Here, the maximal contributions Ch u ~ 15 TeV 2 and ~ —15 TeV 2 . 
Frame c) compares results from the two models using a common scale. Here, it is clearly 
seen that the mSUGRA model is enormously fine-tuned compared to the RNS2 benchmark 
point. 



> 0.01 
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mSUGRA 

H„ 
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H,, 



3 - 



mSUGRA 
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Figure 3: Signed contributions to M§/2 from terms in the EWSB minimization condition Eq. 1.1 
from a) the RNS2 benchmark point defined in the text and 6) the corresponding mSUGRA model 
as RNS2 with fi and tua as outputs rather than inputs. In frame c), the results for both models 
are plotted on a common scale. 

Our goal in this paper is to provide a rather complete characterization of radiative 
natural SUSY. This should provide a comprehensive picture as to where in model parameter 
space we can find 1) ~ 125 GeV along with 2) low EWFT Aew 30 while at the same 
time 3) respecting LHC constraints on sparticle masses. With this goal in mind, in Sec. 2 
we show parameter space regions leading to low Aew from scans over the 2-parameter non- 
universal Higgs model NUHM2 which allow for radiative natural SUSY. In Sec. 3 we extend 
the results to include the split generation non-universal Higgs model NUHM3, wherein 
high scale third generation scalar masses mo (3) need not equal first/second generation 
scalar masses mo(l,2). While the former implementation allows for fewer parameters, the 
additional freedom in the NUHM3 model allows for a more robust decoupling solution 
to the SUSY flavor and CP problems because heavier multi-TeV first /second generation 
sfermion masses are then possible. In Sec. 4, we show that constraints from £>-physics - 
especially BF{b — > 57) are much more easily respected in RNS than in generic NS models. 
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In Sec. 5, we discuss prospects for detecting RNS at the LHC. We also show the new RNS 
mo vs. m\/2 parameter plane, which offers a template for future searches for RNS at the 
LHC. Searches for RNS at the ILC are discussed in Sec. 6 while direct and indirect detection 
of higgsino-like WIMPs is discussed in Sec. 7. In an Appendix, we present formulae needed 
for implementation of our measure of electroweak fine-tuning Aew- 

2. Radiative natural SUSY from the NUHM models 

The direct supersymmetrization of the Standard Model - augmented by weak scale soft 
supersymmetry breaking terms - leads to the Minimal Supersymmetric Standard Model 
(MSSM). Since the mass scale of the MSSM is stable to radiative corrections even when 
the MSSM is embedded into a high scale framework, it is tempting to speculate that the 
MSSM arises as the low energy limit of an underlying SUSY grand unified theory with 
a unification scale Mqtjt — 2 x 10 16 GeV. Indeed, the MSSM (possibly with additional 
gauge singlets and/or additional complete SU(5) multiplets) receives some indirect support 
from experiment in that 1) the measured weak scale gauge couplings nearly unify at Mqtjt 
under MSSM RG evolution, 2) radiative corrections due to the large top quark Yukawa 
coupling - consistent with mt ~ 173 GeV - dynamically breaks electroweak symmetry, and 
3) a light SM-like Higgs boson has been discovered to be lying squarely within the narrow 
mass window predicted by the MSSM. 

Motivated by these successes, the interesting question arises as to whether a natural 
SUSY sparticle mass spectrum, i.e. one with a modest value of Aew> can be consistently 
generated from a model with parameters defined at the high scale Q = Mgtjt- Natural- 
ness requires ~ Mz -^/Aew/2, whilst the recently measured[41] value of the branching 
fraction BF(B S — > qualitatively agrees with the predicted SM value, which in turn 

requires the CP odd boson A to be relatively heavy. We are thus led to adopt the 2- 
parameter non- universal Higgs model (NUHM2)[32], wherein weak scale values of ji and 
rriA may be used as inputs in lieu of GUT scale values of m 2 Hu and m 2 Hd . 2 In Sec. 2.1, 
for simplicity we take a common GUT scale mass parameter mo for all the matter scalars. 
Motivated by Grand Unification, we also assume a common GUT scale gaugino mass pa- 
rameter. Later in Sec. 3 we also explore the possibility of split first/second versus third 
generation matter scalars where we allow the third generation GUT scale mass parameter 
mo(3) to differ from mo(l,2) for the first/second generation scalars. Universality within 
each generation is well-motivated by 50(10) GUT symmetry, since all matter multiplets of 
a single generation belong to a 16-dimensional spinor rep of SO(W). We can also envisage 
some degree of non-universality between mo(l) and m-o(2) as long as both lie in the tens 
of TeV regime: such a scenario invokes a partial decoupling-partial degeneracy solution 
to the SUSY flavor and CP problems (for constraints from FCNC processes [42], see e.g. 
Ref. [43]). For convenience, we will take mo(l) = mo(2). 

2 Since the Higgs fields belong to different multiplets from matter fields, it is easy to envisage models 
with independent SUSY breaking mass parameters for Higgs and matter scalars. 
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2.1 RNS from the NUHM2 model 



The NUHM2 model is defined by the NUHM2 parameter set (1.8) introduced earlier. We 
take m t = 173.2 GeV throughout this paper. For our calculations, we use the Isajet 
7.83 [44] SUSY spectrum generator Isasugra[45]. Isasugra begins the calculation of the 
sparticle mass spectrum with input DR gauge couplings and f T Yukawa couplings 
at the scale Q = Mz (ft running begins at Q = mt) and evolves the 6 couplings up 
in energy to scale Q = Mgut (defined as the value Q where gi = 52) using two- loop 
RGEs. We do not enforce the exact unification condition g% = g\ = 52 at MquTj since a 
few percent deviation from unification can be attributed to unknown GUT-scale threshold 
corrections [46]. Next, we use the SSB boundary conditions at Q = Mgut and evolve 
the set of 26 coupled two-loop MSSM RGEs [47, 48] back down in scale to Q = M z . 
Full twodoop MSSM RGEs are used for soft term evolution, and the gauge and Yukawa 
coupling evolution includes threshold effects in the onedoop beta-functions, so the gauge 
and Yukawa couplings transition smoothly from the MSSM to SM effective theories as 
different mass thresholds are passed. In Isasugra, the values of SSB terms which mix 
are frozen out at the scale Q = MgusY = ^ m t L m i R i while non-mixing SSB terms are 
frozen out at their own mass scale [45]. The scalar potential is minimized using the RG- 
improved one-loop MSSM effective potential evaluated at an optimized scale Q = Msusy 
to account for leading two- loop effects [49]. Once the tree- level sparticle mass spectrum 
is obtained, one-loop radiative corrections are calculated for all sparticle and Higgs boson 
masses, including complete one-loop weak scale threshold corrections for the top, bottom 
and tau masses at scale Q = MgusY [50]. Since Yukawa couplings are modified by the 
threshold corrections, the solution must be obtained iteratively, with successive up-down 
running until a convergence at the required level is found. Since Isasugra uses a "tower of 
effective theories" approach to RG evolution, we expect a more accurate evaluation of the 
sparticle mass spectrum for models with split spectra than with programs such as SuSpect, 
SoftSUSY or Spheno, which make an all-at-once transition from the MSSM to SM effective 
theories. 

Our goal in this section is to find parameter ranges of the NUHM2 model which 
satisfy LHC sparticle and Higgs boson mass constraints while maintaining a low level of 
EWFT. We will also calculate the allowed mass range for various sparticles in low fine- 
tuned/phenomenologically viable parameter space. Toward this end, we search for regions 
of the NUHM2 parameter space with Aew ^5 30, where fine-tuning is better than about 
3%. We will also require that our calculated light Higgs scalar mass lies within the range 
m/i = 125 ± 2 GeV to allow for an estimated uncertainty in our calculation of m,h- We will 
also require that the parameters and m^ 2 respect the recent LHC limits on squark and 
gluino masses obtained within the mSUGRA model [4, 5]. 

We search for radiative Natural SUSY solutions by first performing a random scan 
over the following NUHM2 parameter ranges: 

m : 0-20 TeV, 
m 1/2 : 0.3-2 TeV, 
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-3 < Ao/mo < 3, 
H : 0.1 - 1.5 TeV, 



(2.1) 



0.15 - 1.5 TeV 



tan /3 



3 - 60. 



We require of our solutions that: 

• electroweak symmetry be radiatively broken (REWSB), 

• the neutralino Z\ is the lightest MSSM particle, 

• the light chargino mass obeys the model independent LEP2 limit, m^ > 103.5 GeV[51], 

• LHC search bounds on m g and nig are respected, 

• m h = 125 ± 2 GeV. 

To begin our investigation of NUHM2 model parameters leading to low Aew > m Fig- 4 
we plot each scan point as a red "+" in frames of Aew versus a) mo, b) m 1 / 2 , c) Aq/tuq, 
d) tan/3, e) /i, and f) tua- Since low Aew solutions are only possible for low values of 
fx, we have performed a separate narrow scan, but this time with /i restricted between 
100-300 GeV. The results of this second scan is shown by the blue crosses in the figure. 

We see from the plots that Aew varies from as low as ~ 10 (A^ = 10% EWFT) 
to over 1000. While the bulk of points shown are fine-tuned with large Aew ^ 100, there 
do exist many solutions with Aew i$ 30, corresponding to better than 3% EWFT. The 
RNS solutions with Aew 30 are obtained for values of too ~ 1 — 8 TeV. In the cases 
where too is as high as 5-10 TeV, the top squark masses are driven to much lower values via 
1) the large top-quark Yukawa coupling ft which suppresses top-squark soft masses during 
RG evolution, 2) large mixing effects which can suppress and yield a large — m^ 2 
splitting, and 3) two-loop RGE suppression of diagonal top squark mass terms arising from 
large first/second generation sfermion masses[52, 53, 54]. If iuq is too large - in this case 
above ~ 10 TeV - then these suppression mechanisms are insufficient to drive m^ to low 
enough values to allow for low EWFT. Thus, the span of points shown in frame a) trends 
upward in Aew as mo increases past about 8 TeV. We also see that for the red pluses in 
frame a) Aew has an upper bound close to about 500 if too < 10 TeV. For still larger values 
of itiq then Aew increases with too- This is because while fi 2 (or equivalently — m\ ) is 
the largest of the quantities in (1.1) for the lower range of too, for very large values of mo 
then E" begins to dominate. The blue crosses from the narrow scan with small \x have a 
different shape from the red broad scan since the upper edge is mostly determined by S^, 
and so increases with mo- 

In frame b) of Fig. 4, we show Aew vs - m i/2- Here, the low values of Aew span a wide 
range of m^ 2 values from 0.3 — 1.5 TeV. Since m~ g ~ (2.5 — 3)m 1/ / 2 , we expect Aew ^5 30 for 
rag values up to about 4 TeV. For the entire parameter space (red pluses) Aew is roughly 
evenly distributed with respect to the gaugino mass parameter. In frame c), we show 
Aew vs - Ao/mQ. We see a clear trend for low values of EWFT when |^4o/ m o| ~ 1-5 — 2. The 
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Figure 4: The dependence of Aew on various NUHM2 parameters from a scan (2.1) over parameter 
space (red pluses) and for the dedicated scan with 100 GeV< ^ <300 GeV (blue crosses). The line 
at A EW = 30 is to guide the eye. 

reason is that the hole at low magnitudes of Aq / mo and small values of Aew occurs because 
of the Higgs mass constraint. Large magnitudes of GUT scale Aq lead to correspondingly 
large weak scale At parameters, which, in turn, provide large mixing in the top-squark 
sector. This leads to low EWFT and also heightened values of rrih ~ 125 GeV. Frame d) 
shows Aew vs - tan/3. We see a slight preference for low tan/3 ~ 10 — 20 but otherwise 
no structure to speak of. Frame e) shows Aew versus the weak scale value of /i. The 
parabolic lower edge of the span of points reflects the upper bound on [i necessary for low 
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EWFT. From the plot, bounds on \x can be conveniently read off: for instance, requiring 
Aew i$ 30 then requires fi < 350 GeV. Of course, models with low fi ~ 100 GeV but 
multi-TeV top squarks can still be very fine-tuned if the dominant contributions to Aew 
arise from S^(tj). In frame f), we plot Aew vs. vtia- We see that low Aew can be found 
over the entire range of tjia ~ 0.15 — 1.5 TeV, so this parameter is not so relevant towards 
achieving low EWFT. 

Next, to gain a sense of the sparticle mass ranges expected from RNS, we plot Aew 
versus selected sparticle masses. First, since tuq ~ 2 — 8 TeV for Aew ^5 30, we expect 
first and second generation squark and slepton masses also within this range (which is for 
the most part inaccessible LHC SUSY searches). Next, in frame a) of Fig. 5, we show 
Aew vs - m g- We find that requiring Aew ^5 30 requires rrig ~ 1 — 4 TeV. The lower 
portion of this range should be accessible to LHC14 searches, while the upper part lies 
beyond any LHC luminosity upgrade[55]. 




2 4 6 8 10 12 14 16 18 2 4 6 8 1 1 2 14 1 6 1 8 

m- u (TeV) m s , (TeV) 



Figure 5: The value of Aew versus gluino and third generation squark masses from a scan over 
NUHM2 parameter space. As in Fig. 4, the red pluses denote the distributions from the complete 
scan, whereas the blue crosses depict the results for the dedicated low [i scan. The line at Aew = 30 
is to guide the eye. 

In frame b), we show Aew versus the lighter top squark mass mj- . Here, we see that 
Aew ^5 30 allows mr ~ 0.5 — 2.5 TeV range. This is well above the range expected in 
generic NS models[9, 10], where rn^ 12 has been advocated to lie below about 600 GeV. 
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In frame c), we show Aew vs. mj 2 . Here, we find that mf 2 can range up to ~ 6 TeV 
for Aew ^5 30. Such high values of m t - are helpful to increase radiative corrections to 
the light Higgs mass into the 125 GeV range. However, such heavy top squarks lie 
far beyond any conceivable LHC reach. In frame d), we show Aew vs - m b 1 - Here we see 

~ 0.8 — 6 TeV, which again allows for far heavier bottom squarks than previous NS 
models, where 2 and mg all were suggested to be < 600 GeV. 

In Fig. 6a), we show Aew vs - m S^ x - For RNS models, — m z 1 2 ~ M> ^ s mass 
is roughly equal to that of the two lighter neutralinos. Since W\ is mainly higgsino-like 
near the lower edge of the envelope of points, the distribution follows a similar pattern as 
for the Aew vs. fi plot in Fig. 4. We see that for Aew ^ 20, < 250 GeV. Thus, a 

linear collider operating with ^fs > 2m^ will directly probe the lowest (and hence most 
lucrative!) values of Aew if the relatively soft visible daughters of the chargino can be 
distinguished over two-photon backgrounds [5 6]. In this sense, it has been emphasized that 
for models of natural SUSY, a linear e + e~ collider would be a higgsino factory in addition 
to a Higgs factory[14, 11, 20]! In frame b), we show Aew vs - m w 2 ' ^ n ^ e RNS model, the 
W2 is nearly pure wino-like and its mass can range between ~ 0.3 — 1.2 TeV for Aew ^ 30. 
Since RNS as presented here includes gaugino mass unification, then typically Zip are 
higgsino-like, Z3 is bino-like and Z^ is wino-like. Since the SU(2) gauge coupling g is 
rather large, we expect significant rates for W^Z\ production at LHC, at least for the 
lower portion of the range of ' ^ n f rame c ) 5 we show the — m z 1 mas s difference in 
RNS versus Aew- For most points with Aew % 30, we find that — < 10 — 20 GeV. 
Some points with Aew ~ 30 — 40 have a mass difference as large as 100 GeV; these points 
arise from sampling the lower portion of the mi 12 range, which gives rise to gaugino masses 
comparable in magnitude to so that the lighter electroweakinos are actually gaugino- 
higgsino mixtures. For the more likely small mass gap case, the lighter neutralinos are 
dominantly higgsino-like and decay via Z2 —> Ziff (here / denotes SM-fermions) through 
the virtual Z. Then decays into opposite-sign same-flavor (OS/SF) isolated dileptons 
should occur at ~ 3% for each charged lepton species. The presence of low invariant mass 
OS/SF isolated dileptons from boosted Z2 produced in gluino or gaugino cascade decay 
events could then be a distinctive signature of RNS at the LHC. For NUHM2 models with 
larger values of Aew> falling outside the RNS low EWFT requirement, can be greater 
than tti^ + Mz or + irih so that two body decays of Z2 are then allowed. Finally, in 
frame f), we show Aew vs - m h- Here, we see the lower ~ 123 — 124 GeV values are 
just slightly preferred by EWFT over the higher range, although values of rrih as high as 
~ 126.5 GeV occur for A EW = 30. 

While our methodology allows one to find a low value of /i 2 for any value of mo and 
mi/2, this by itself does not guarantee a small value of Aew- in addition, the GUT scale 
value of rrvjj has to be adjusted very precisely to obtain low EWFT, which could be viewed 
as a different sort of fine-tuning: that only a very narrow range of m 2 Hu {MQui) values will 
yield —m 2 Hu ~ M§ at the weak scale. 3 This can be seen from Table 1 where we plot the 

3 From the perspective introduced in Sec. 1, we would look for an underlying model where is thus 

determined. 
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Figure 6: The value of Aew versus electroweak -ino and Higgs boson masses from a scan over 
NUHM2 parameter space. As in Fig. 4, the red pluses denote the distributions from the complete 
scan, whereas the blue crosses depict the results for the dedicated low /i scan. The line at Aew = 30 
is to guide the eye. 

value of Tn/jj (Mqxjt) which is needed to generate small \x solutions for two different cases 
of NUHM2 model parameters. Optimistically speaking, we would view this as essentially 
determining the GUT scale value of m 2 Hu /m^ to be very nearly 1.65 (Case A) or 1.71 (Case 
B). It is gratifying to see that the GUT scale values of all scalar mass parameters have 
no hierarchy as expected in models of gravity-mediated SUSY breaking where all scalar 
masses might be expected to be comparable at the high scale. 



Case A 




Case B 




m 2 Hu (M GU t) 




Aew 


m 2 Hu (M GU t) 




Aew 


1.03 x 10 7 


150 


9.04 


2.73 x 10 7 


150 


15.4 


1.02 x 10 7 


250 


18.8 


2.72 x 10 7 


250 


24.1 


1.00 x 10 7 


400 


42.4 


2.70 x 10 7 


400 


49.5 



Table 1: An illustration of the sensitivity of the EWFT fine-tuning measure Aew to m 2 H (Mqut)- 
For case A the NUHM2 parameters are m = 2.5 TeV, m 1/2 = 400 GeV, tan/3 = 10, m A = 1 TeV, 
while for case B we have, m = 4 TeV, m 1 / 2 = 1 TeV, tan/3 = 15, m A = 2 TeV. For both cases, 
we take Aq = — 1.6mo. The numbers in the Table are in GeV units. 
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To be more general, we show in Fig. 7 a scatter plot of Aew versus the GUT scale ratio 
m 2 H /mo from our scan over NUHM2 models. We find that for points with Aew % 30, 
then m Hu (M G u T ) ~ (1 - 2)m . 




Figure 7: The value of Aew versus m# u /mo(MGUT) from the scan over the NUHM2 parameter 
space. As before, the red pluses arc for the scan over the entire range of fi while the blue crosses 
are for the dedicated scan with /i limited to the 1 00-300 GeV range. The line at A^w = 30 is to 
guide the eye. 

2.2 RNS from the NUHM1 model? 

Up to this point we have focused on RNS from the NUHM2 model. However, it is of interest 
to see if low EWFT is also possible within the NUHM1 [32] framework in which H u and 
Hd have equal GUT scale mass parameters; i.e. m^(MGUT) = ^^(-^gut) = fr^. 
As mentioned above, m? is then adjusted to that m^(MsusY) ~ M|. For brevity, we 
confine our investigation to model lines where we fix mo, m^j and Aq to be the same as 
for the RNS2 model point, but where we vary m^, the common GUT scale Higgs mass 
parameter and tan/3. In Fig. 8 we show the value of a) fj,, b) uia, and c) Aew versus 
m^ for tan/3 = 8.85 (the RNS2 value), 25, 40 and 50. We see that the various curves 
in frame a) are quite close (except at very large values where they dive down). This is 
essentially because the top Yukawa coupling that dominantly affects vrr H (remember that 
mjj (-Mstjsy) determines fi) hardly varies with tan/3; the small differences arise from the 
(subdominant) effects of bottom Yukawa couplings. In contrast, the tua values in frame 
b) reduce considerably as tan/3 increases. We can understand this if we remember that 
the bottom Yukawa coupling- which increases with tan/3- drives m'j{ d to low values thus 
reducing m 2 A ~ rrijj + [J? for larger tan ft values. Turning to frame c) we see that for this 
model line with tan/3 = 8.55 (uppermost curve), Aew reduces with increasing m^ as in 
the fi curve in frame a) until the kink at which it starts increasing. We have checked that 
the kink occurs when fj? becomes so low that the frijj term becomes larger than all other 
terms in (1.1). For larger values of tan/3, the m 2 H contribution is suppressed, resulting in 
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smaller values of Aew- However, in none of the cases shown does Aew drop below ~ 80. 
It may be a useful exercise to scan the NUHM1 parameter space to see just how small the 
EWFT can be when all LHC constraints are satisfied. 





Figure 8: Plot of a) fi, b) itia and c) Aew in the 1-parameter NUHM1 model versus for RNS2 
model parameters mo = 7025 GeV, TOx/2 = 568.3 GeV, A = —11426.6 GeV and for several values 
of tan/3. In all the frames the order of the lines is that of increasing tan/3, with tan/3 = 8.55 on 
the top and tan /3 — 50 on the bottom. 



3. RNS from the NUHM3 (split generation) model 

In this Section, we investigate if any advantage can be gained for RNS models if we allow for 
a splitting between scalars of the third generation and those of the first/second generations. 
We adopt the parameter set 

m (l,2), m (3), m 1/2 , A , tan/3, fj,, m A (NUHM3) (3.1) 

where mo (3) is the GUT scale third generation soft SUSY breaking mass parameter and 
mo(l,2) is the corresponding (common) parameter for the first/second generation. 
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We search again for RNS solutions from the split generation parameter space by im- 
plementing a random scan over the parameters: 



m (3) : 0-20 TeV, 
m (l,2) : m Q (3) -30 TeV, 



»™l/2 



0.3 - 2 TeV 



3 < Aq/ttio < 3, 



(3.2) 



0.1 - 1.5 TeV 



m A 



0.15 - 1.5 TeV 



tan/3 



3-60. 



We implement the same LHC sparticle mass and rrih = 125 ± 2 GeV constraints as before. 

In Fig. 9, we show Aew versus mo (3) and also versus mo(l,2). The results for Aew 
versus other parameters are very similar to Fig. 4 so we do not repeat them here. From 
Fig. 9a), we see that RNS solutions with Aew % 30 can be found for mo (3) values ranging 
between 1-8 TeV, similar to results found in Fig. 4 for the NUHM2 model. It is interesting 
to note that the smallest values of Aew i n the figure are no smaller than for the NUHM2 
model. The gap at small values of mo (3) is an artifact of the upper limit on m^ m our 
scan: for small values of mo (3) the lighter t-squark is often driven to tachyonic masses 
by two-loop contributions of heavy first/second generation squarks. We have checked that 
with larger values of m^ 2 in the scan, solutions fill in the entire gap. Again, even though 
the GUT scale value of mo (3) is in the multi-TeV regime, the t-i and especially t\ physical 
masses are considerably lower - in the few TeV regime - due to radiative effects from RGE 
running and also large mixing. 

The key advantage of the NUHM3 model is seen in Fig. 96), where we plot Aew 
versus mo(l,2). In this case, we see that GUT scale first/second generation scalar masses 
can easily range between 1 — 30 TeV while still maintaining low Aew- The solutions 
with mo(l, 2) in the multi-TeV region will also produce first/second generation squark and 
slepton masses which are comparable to mo (1,2). This allows for a much more robust 
solution to the SUSY flavor/CP problems. It also ameliorates the cosmological gravitino 
problem if m.3/2 ~ mo(l, 2) as is expected in simple models of gravity- mediation. 

We do not show plots of Aew versus sparticle masses since these are very similar to 
results shown in Fig. 5 and Fig. 6 except for the fact that NUHM3 scans allow for much 
heavier first/second generation squark and slepton masses in the 10-30 TeV range, whereas 
in the NUHM2 model the squarks and sleptons are typically constrained to be below 8 TeV 
due to the imposed relation mo(3) = mo(l,2). 



4. Rare B decay constraints on RNS 
4.1 BF(6 -»• si) 

The combination of several measurements of the b — > decay rate finds that BF(b — > 
s 7 ) = (3.55 ± 0.26) x 10 -4 [57]. This is slightly higher than the SM prediction [58] of 
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Figure 9: The value of Aew versus mo(3) and mo(l, 2) from a scan over NUHM3 model with split 
first/second and third generations. As in Fig. 4, the red pluses denote the distributions from the 
complete scan, whereas the blue crosses depict the results for the dedicated low fi scan. The line 
at Aew = 30 is to guide the eye. 

BF SM (b -)• S7) = (3.15 ± 0.23) x 1CT 4 . SUSY contributions to the b ->■ S7 decay rate 
come mainly from chargino-stop loops and the W-charged Higgs loops, and so are large 
when these particles are light and when tan/3 is large[59]. Thus, in generic natural SUSY 
where g < 500 GeV, one finds generally large deviations from the SM value for 
BF(b 57), as shown in Ref. [11]. In contrast, in radiative natural SUSY where third 
generation squarks are in the TeV range, SUSY contributions to BF(b — > S7) are more 
suppressed. The situation is shown in Fig. 10a) along with the measured central value 
(green solid line) and errors. The red points all have Aew < 30 and qualify as RNS points. 
We see the bulk of RNS points are consistent with the measured BF(b — > sj), although 
there are outliers. 




BF(B s -> ji \i) 



Figure 10: The values of Aew versus a) BF(b —> sj) and b) BF(B S ). The vertical lines 

represent the experimental measurements with uncertainties. 

4.2 B s fi+i^r 

Recently, the LHCb collaboration has discovered an excess over the background for the 
decays /i + /i"[41]! They find a branching fraction of BF(B S fi + fi~) = 3.2i};^xl0~ 9 
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in accord with the SM prediction of (3.2 ± 0.2) x 10 -9 [60]. In supersymmetric models, 
this flavor-changing decay occurs through pseudoscalar Higgs A exchange[61], and the 
contribution to the branching fraction from SUSY is proportional to ( tan ^) . "We show the 

m A 

value of BF(B S — > from RNS in Fig. 106). The decay is most constraining at large 

tan/3 ~ 50 as occurs in Yukawa-unified models[62] and low niA- In the case of RNS with 
lower tan/3 and heavier m^, the constraint is less important. The bulk of the RNS points 
in Fig. 106) fall well within the newly measured error bands although there are some outlier 
red points, mainly at larger values of the branching fraction. 

4.3 (g - 2)„ 

In addition, the well-known (g — 2)^ anomaly has been reported as a roughly 3a deviation 
from the SM value: Aa M = (28.7 ± 8.0) x 10 _10 [63]. In RNS, since the [L\^ and masses 
are expected to be in the multi-TeV range, only a tiny non-standard contribution to the 
{g — 2)^ anomaly is expected, and alternative explanations for this anomaly would have to 
be sought. 



5. RNS at LHC 



Here we list a few of the possibilities for an LHC search for radiative natural SUSY. A 
thorough study of signal and background simulations will be presented in an upcoming 
study [64]. 

The hallmark feature of radiative natural SUSY models is the presence of light higgsino 
states W\ and Z\p. with masses ~ \p\ ~ 100 — 300 GeV, and usually a small mass gap 
— frig and — of order 10-30 GeV with a possible exception of low m 1 / 2 and 
larger /i where there can be substantial gaugino-higgsino mixing. 

One possibility for RNS at LHC is to search for clean trilepton events from W\Z 2 
production followed by W\ — > Z\£vi and Z2 — > Z\i + £~ decays where i = e or /i. This 
signal has been investigated in Ref. [14]. There, the Pt(£) values were typically found to 
be quite low in the 5-15 GeV range making detection difficult. The small mass difference 
between the parent and daughter neutralino will also mean that the invariant mass of the 
opposite-sign same-flavor dilepton pair will be small, making it more challenging to separate 
it from SM origins. Nevertheless, this reaction certainly motivates an LHC search for clean 
trilepton states with very soft lepton values, as low as is experimentally feasible. 

While first/second generation squarks are expected to be in the multi-TeV range, 
the value of m g is expected to be ~ 1 — 5 TeV. The lower portion of this mass range 
m g ~ 1 — 2 TeV should be accessible to LHC searches for gluino pair production pp — > gg. 
For RNS models, since nig, then gluino three-body decays to third generation 

particles typically dominate: g — > tbW\ or tiZi. Thus, we would expect gg events to contain 
up to four 6-jets, and 2-4 reconstructable top quarks. A small fraction of events would 
contain Z2 — > Z\l + t~ where m(£ + £~) is bounded by — ~ 10 — 20 GeV. Normally 
the leptons from Z% decay would be rather soft, but in the case of large boosts from the 
gluino cascade decay, the opposite-sign/same-flavor pair would be highly collimated in 
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opening angle. We expect LHC14 with 100 fb to be able to probe rrig ~ 1 — 2 TeV via g 
cascade decays using analyses similar to those used for gluino searches in mSUGRA when 
mo is very large. [55]. 

A novel search for RNS at LHC is to look for pair production of the heavier gaugino 
states W2 and Z3 and Z4. Wino pair production occurs via the large SU(2) gauge couplings 
and leads to large rates for W2Z4 and W2W2 processes. The decays W2 — > WZ\^ZW\ 
and Z4 — > WW\, ZZ\ t 2 occur with significant branching fractions and yield a variety of 
diboson final states that include spectacular W^W^ 1 and WZ [65] plus E™ 1SS events with 
soft debris from the decays of higgsinos. Events with light Higgs bosons instead of gauge 
bosons in the final state are also possible [66]. 

As with the mSUGRA model, a wide range of RNS signatures for LHC can be found 
by exploring the mo vs. 7771 / 2 plane. This plane will look quite different from the mSUGRA 
case since now we will require small /j ~ 100 — 200 GeV in accord with EWFT and also 
Aq ~ — 1.67770 in accord with rrih = 125 GeV and low EWFT. We show the plane in Fig. 11 
for u, = 150 GeV with Aq = — 1.67770, tan/3 = 10 and 777,4 — 1 TeV. Here, we plot contours 
777/j = 123 and 125 GeV and also contours of Aew = 6, 10, 15 and 50. Almost the entire 
plane has low Aew < 50, with 5.5 < Aew 10 in the lower left portion. In addition, the 
right-hand portion of the plane has m/i > 123 — 125 GeV. The purple-shaded region marked 
LEP2 has < 103.5 GeV in violation of LEP2 limits on chargino pair production. We 
also show recent LHC constraints from gluino/squark searches within the mSUGRA model 
as the black contour. 4 We extrapolate those constraints to much higher mo values via the 
dashed black contour. Over the entire plane, 777^ ~ m z 1 ~ /j = 150 GeV so there would 
always be light higgsino pair production at LHC. The region with 777i/ 2 ^$ 0.6 TeV yields 
777<j < 2 TeV and should be accessible to future gluino pair production searches. Signals 
from wino pair production may also be observable at the LHC, and perhaps even at LHC8, 
if the heavier chargino is sufficiently light. 

6. RNS at ILC 

Since the main feature of RNS models is the presence of light higgsinos W\, Z\ and Z2, we 
expect excellent prospects for testing RNS at a linear e + e~ collider. Pair production for 
charged higgsinos via the reaction e + e~ — > W^W^ would yield soft but observable decay 
products from W\ — > Z\f f decay, where / and /' are SM fermions. These decay products 
should be easily detectable in the clean environment of e + e~ colliders, and moreover should 
be acollinear in the transverse plane as opposed to two-photon backgrounds 77 — > ff where 
visible decay products tend to come out back-to-back. The entire mo vs. plane shown 
in Fig. 11 will be accessible to an ILC with y/s > 2m^ ~ 2|/x|. 

Moreover, the cross section for the reaction e + e~ — > Z1Z2 should also be large and 
provide corroborative evidence. These cross sections will have a distinctive shape versus 

4 Strictly speaking, these are the constraints obtained in the mSUGRA model for Aq — from the 
non-observation of signals from gluino and first generation squark production. Since the masses of these 
sparticles depend mostly on mo and m 1( i 2 , it is reasonable to suppose these also apply to the plane in 
Fig. 11. 
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NUHM2: tan$=10, A Q =-1.6m , yx =150 GeV, m t =173.2 GeV 
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Figure 11: Contours of Aew ( re d curves) and (purple curves) in the Too versus m 1 / 2 plane of 
RNS model with A = -1.6m , tan/3 = 10, n = 150 GeV and m A = 1 TeV. 

beam polarization as shown in Ref. [14] which should be indicative of higgsino pair pro- 
duction. We note here that since Ecu ^ ^ rn Wi ~ 2 1 /^t | , and Aew ~ M then an 
e + e~ collider with a CM energy Eqm directly probes 



so that even a low energy e + e~ linear collider would probe the most lucrative regions 
of RNS parameter space (that portion with lowest Aew) an d as Ecm increases, would 
discover natural SUSY or increasingly exclude it. 

7. Search for higgsino-like WIMPs from RNS 

One of the distinctive features of natural SUSY models is that the lightest MSSM particle 
is a higgsino-like neutralino Z\. If i?-parity is conserved, then the Z\ may make up all 
or at least a portion of the dark matter in the universe. Higgsinos with mass > 
mw, Mz have high annihilation rates into vector boson pairs. Thus, if they are present in 
thermal equilibrium in the early universe, then the higgsino relic density may be computed 
approximately as 



z i p c /h 2 \Trg*J mpi (crv) 
where sq is the entropy density of the universe at the present time, p c is the critical closure 
density, h is the scaled Hubble constant, g* is the number of relativistic degrees of freedom 



Aew ~ E 2 CM /{2M 2 Z ) 



(6.1) 




(7.1) 
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at freeze-out, Xf ~ 25 is the scaled freeze-out temperature, mpi is the Planck mass and 
(av) is the thermally averaged neutralino annihilation cross section times relative velocity. 
Higgsino-like WIMPs couple with gauge strength to vector bosons so that (av) is large and 
the relic density is suppressed. 

We evaluate the relic density of higgsinos using Isatools[67] from our scan over NUHM2 
parameters as in Sec. 2 and show Aew vs - * n f rame a ) °f Fig. 12. Points with 

Aew < 30 indicative of RNS are shown in red while the more highly fine-tuned points 
are in blue. The vertical green line shows the WMAP-measured value of the dark matter 
density. We see that the majority of points with Aew < 30 have ^j^/i 2 ~ 0.005 — 0.05, 
i.e. well-below the measured abundance. Several points have ^■^ 1 h 2 > 0.12; these points 
arise from cases where [i ~ M\ where the neutralino is of mixed bino-higgsino variety. 

There exist a variety of non-standard cosmologies with features which make them more 
attractive than the standard WIMP-only dark matter scenario. For instance, in stringy 
models with moduli fields at the 10-100 TeV scale, the moduli may decay after BBN into SM 
particles, thus diluting all relics present. Alternatively, if moduli decay to SUSY particles 
which cascade into the LSP, then the neutralino abundance may be enhanced[68]. 

Another possibility arises from SUSY models where the strong CP problem is solved by 
the Peccei-Quinn mechanism[69] with its concomitant axion a. In the SUSY case, the axion 
superfield a also contains an i?-even spin-0 saxion s and an i?-odd spin-^ axino a. In the case 
where Z\ is LSP, then the dark matter would be a mixture of two particles: the axion and 
the neutralino. Axinos which are produced thermally at high re-heat temperature Tr in the 
early universe would cascade decay to neutralinos at a decay temperature Tjj < TV , causing 
a neutralino re- annihilation which provides a much higher abundance of WIMPs than 
expected in a WIMP-only picture. In addition, saxions can be produced both thermally 
and via coherent oscillations, and may decay to both SUSY and SM particles: the former 
case enhances the neutralino abundance while the latter case dilutes any relics present at 
the time of decay. Calculations of the neutralino abundance in the PQ-augmented MSSM 
depend on the various PQ parameters along with Tr and the SUSY particle spectrum and 
have been presented in Ref. [70, 71, 72]. In the case of models with a standard under- 
abundance of neutralinos, the neutralino abundance is almost always enhanced beyond its 
standard value Vl-^h?. If this scenario is applied to the case of RNS, then we may most 
likely expect an enhanced higssino-like WIMP abundance beyond its standard value. In 
this scenario, axions will also be produced via coherent oscillations at temperature around 
the QCD phase transition. Thus, the higgsinos could make up either a small or a large 
fraction of the relic dark matter, with axions comprising the remainder. The important 
point here is that it is very difficult to suppress the higgsino abundance below its standard 
thermal value which is shown in Fig. 12a). Thus, we would expect relic higgsinos to be 
present in the universe today, but with an abundance which is suppressed by between 1 — 15 
from the measured value. This opens up the opportunity to detect relic higgsinos, albeit 
while these would only constitute a fraction of the measured dark matter abundance. At 
the same time, there is also the possibility to detect relic axions. 

With a view towards detecting relic higgsinos from RNS, we show in Fig. 126) the value 
of Aew versus the spin-independent neutralino-proton scattering cross section a SI (Zip) in 



- 24 - 




Figure 12: The value of Aew versus neutralino relic density and direct and indirect WIMP 
detection rates. Vertical green lines denote upper experimental limits obtained assuming that the 
WIMP saturates the observed density of dark matter. The predictions in the last three frames need 
to be re-scaled by a factor VL^Ji 2 if the neutralinos make up only part of the dark matter. 



pb from Isatools[73]. The red points with Aew < 30 occur witha SI (Z lP ) ~ Kr 9 -1(T 7 pb. 
For comparison, we show via the green vertical line the Xe-100 limit from 225 live days[74] 
for ~ 150 GeV. Naively, many of the RNS points would be excluded if higgsinos 
comprised the entire dark matter density. However, in the mixed axion/higgsino dark 
matter scenario, the expected local abundance of WIMPs can be scaled down by factors of 
1 — 15 typically. Even with this rescaling of the expected local abundance, we still expect 
relic higgsinos to be within detection range of near- future WIMP detectors. 

In Fig. 12c), we show the spin-dependent neutralino-proton cross section a SD (Z\p) in 
pb. The bulk of RNS points with Aew < 30 populate the region with a SD {Z\p) ~ 2 x 
10~ 5 — 10~ 3 pb. The IceCube neutrino detector at the South Pole is sensitive to detection of 
neutrinos arising from higgsino annihilation in the core of the sun. The expected detection 
rate depends on the sun's ability to sweep up neutralinos via Z\p collisions which depends 
mainly on a SD (Z\p). For reference, we also show the current IceCube WW limit [75] at 
a SD {Z\p) ~ 3 x 10 -4 pb. This limit depends on the assumption that WIMPs comprise the 
entire DM abundance, and would need to be rescaled for a mixed axion/higgsino cosmology. 

Fig. 12c?) shows the thermally averaged neutralino annihilation cross section times 
relative velocity, evaluated as v — > 0. This quantity enters linearly into indirect searches 
for neutralino annihilation in the cosmos into 7s or e + , p or D. For the case of RNS, the 
bulk of points with Aew < 30 inhabit the region around (av)\ v ^o ~ 10~ 25 cm 3 /sec. The 
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vertical green line shows the upper limit on the annihilation cross section times velocity 
for very non-relativistic dark matter in dwarf spheroidal satellite galaxies of the Milky 
Way annihilating to W boson pairs obtained by the Fermi collaboration[76], assuming a 
~ 150 GeV WIMP. Models with a larger annihilation cross section would have led to a flux 
of gamma rays not detected by the experiment, assuming a Navarro- Frenk- White profile [77] 
for each dwarf galaxy in the analysis. We see that the Fermi bound might exclude the bulk 
of points assuming higgsinos saturate the DM density. This bound changes rather slowly 
with the WIMP mass, being just a factor of 2 weaker for a WIMP mass of 300 GeV. Further 
searches and improvements by the Fermi-LAT Collaboration and/or the impending AMS 
results should provide more stringent probes of the RNS model. 

8. Summary and conclusions 

Models of natural supersymmetry reconcile the lack of a SUSY signal at LHC with the 
principle of electroweak naturalness. Natural SUSY models are characterized by light 
higgsinos of mass ~ 100 — 300 GeV, three light third generation squarks with mass less 
than about 500 GeV and gluinos of mass less than about 1.5 TeV. First/second generation 
squarks may be much heavier - in the multi-TeV regime - thus avoiding LHC searches 
and providing at least a partial decoupling solution to the SUSY flavor and CP problems. 
Attractive as they are, generic NS models based on the MSSM are at odds with the recent 
discovery of a light Higgs scalar at ~ 125 GeV which requires TeV-scale top squarks along 
with large top-squark mixing. 

We presented here an improved natural SUSY model dubbed radiative natural SUSY [20] , 
or RNS. RNS is a SUSY model based on the MSSM which may be valid all the way up to 
the GUT scale. Thus, it maintains the desirable features of gauge coupling unification and 
radiative electroweak symmetry breaking while avoiding the introduction of extra possibly 
destabilizing gauge singlets or other forms of exotic matter. The main features of the RNS 
model include 1) a low value of superpotential higgsino mass \fj,\ ~ 100 — 300 GeV, and 2) a 
weak scale value of —vn H ~ M|: both these qualities are required to fulfill electroweak 
naturalness at the tree level. The term vr? H is driven to low values radiatively by the 
same mechanism leading to REWSB and depends on a large top quark Yukawa coupling. 
We proceed further by evaluating EWFT at the 1-loop level. In this case, top squark 
masses enter the computation of Aew an d are also driven radiatively to few- TeV values. 
By allowing for large top-squark mixing (]Aq\ ~ (1 — 2)mo), top-squark contributions to 
EWFT are suppressed at the same time as the light Higgs boson mass is uplifted: thus, 
the model reconciles electroweak fine-tuning with ~ 125 GeV all in the context of the 
MSSM valid up to the GUT scale. 

RNS may be realized in the two-parameter non- universal Higgs models NUHM2. In 
this case, low EWFT with Aew 30 can be attained for model parameters which lead to 
a distinctive mass spectrum: 

• light higgsino-like W\ and Z\ t 2 with mass ~ 100 — 300 GeV, 

• gluinos with mass m,g ~ 1 — 4 TeV, 
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• heavier top squarks than generic NS models: ~ 1 — 2 TeV and ~ 2 — 5 TeV, 

• first/second generation squarks and sleptons with mass m~j ~ 1 — 8 TeV. The 

range can be pushed up to 20-30 TeV if non-universality of generations with 
mo(l,2) > mo(3) is allowed. 

The RNS model with the above spectra also fulfills limits from rare £>-decay measurements, 
which can be an Achilles heel for generic NS models with much lighter third generation 
squarks. 

The RNS model can be tested at LHC for m 1 / 2 in the lower portion of its range 
whereupon gluino pair production and/or gaugino pair production may lead to observable 
signals. Many of the associated SUSY events will contain light higgsinos arising from 
cascade decays with soft decay products which may be observable. The case of OS/SF 
dileptons with mass < 10 — 20 GeV would signal the presence of Z 2 — > Z\i + i~ decay. 

Linear e + e~ colliders would likely provide the definitive test of RNS models since pair 
production of charged higgsinos should be easily observable and the lowest energy machines 
will scrutinize the most lucrative parameter choices with the lowest values of Aew- For 
RNS, an ILC type machine would be a higgsino factory in addition to a Higgs factory. 

In RNS, we also expect the presence of higgsino-like WIMPs which have large rates 
for direct and indirect WIMP detection. Since higgsinos are thermally underproduced, we 
expect them to constitute only a portion of the measured dark matter abundance, with 
perhaps axions comprising the remainder. Detectability via WIMP searches will depend 
on the higgsino fraction of the dark matter. 

The many elegant features presented above impel us to regard RNS as the possible 
new paradigm SUSY model. Its consequences for detection at colliders and at dark matter 
detectors merits a high level of scrutiny. 
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A. Radiative corrections to the Higgs potential minimization conditions 

The Higgs portion of the scalar potential in the MSSM is given by 



Vmggs = U trcc + AV, 



(A.l) 



where the tree level portion for the neutral Higgs sector is given by 



Ut rcc = (m 2 Hu + » 2 )\h° u \ 2 + (m 2 Hd + fi 2 )\h d \ 2 

-BMhS + h.c.) + ltf + ! f)M\*-\hW)* 



(A.2) 
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and the radiative corrections (in the effective potential approximation, and using the DR 
regularization scheme, as appropriate for SUSY models) by 



Ay 



log 



my 
Q 2 



(A.3) 



where the sum over i runs over all fields that couple to Higgs fields, mf are the Higgs field 
dependent mass squared values, and q = c co [C c h a , with c co i = 3 (1) for colored (uncolored) 
particles and c c h a = 2 (1) for charged (neutral) particles and Sj is their spin quantum 
number. 

Minimization of the scalar potential allows one to compute the gauge boson masses in 
terms of the Higgs field vacuum expectation values v u and Vd, and leads to the well-known 
conditions that 



B/iVd 
Bfiv u 



(m 2 Hu + ii 2 



where 



and K,d 



9z( v d- v l)) v u + ^u 

{ m H d + ^ + 9z{ V d - v l)) Vd + ^d: 

dAV 



(A.4) 
(A.5) 

(A.6) 

di)/ v 72 ", 9z — (9 2 +g' 2 )/&- By SU{2) invariance, the scalar potential 



V depends on the scalar fields as [6] V(huh u , h d hd, h u hd + c.c), then we have 



^d 



T%v u + Y,iv d and 



(A.7) 
(A.8) 
(A.9) 



where 



dAV 



d\h u \* 
dAV 



and 



d\h d \< 

S d = dAV 
" d(h u h d + c.c. 

In this case, the minimization conditions may be expressed as 

^ 2 + E3) - (m 2 Hu + E£) tan 2 /3 



M|/2 



tan 2 /3 - 1 



5/i = (m^+^ + E^) + (m^+^ + ES sin/3 cos /3 + E*. 



(A.10) 
(A.ll) 
(A.12) 



(A.13) 
(A.14) 



The advantage of writing the minimization conditions in terms of E" and E^ and T, u d is 

2 Hd and B\i are neatly separated so that terms do not 
appear in Eq. (A.13), and so do not contribute to the fine-tuning calculation. 



that the corrections to m 2 H , 
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The contributions of the various Es can be written as 
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where 



F(m ) = m ( lo & g 2 



with the optimized scale choice Q 2 

For the top squark contributions, we find 
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)/2 + M§ cos 2/3(j — |xw) and x^y = sin 2 6\y. In the denominator 
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of (A.17) and (A.18), the tree level expressions of m 2 ^ should be used. 



For 6-squark contributions, we have 
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where A b = (m- —m~ )/2-M| cos2/3(|-^x H /)- The expressions for £^(fi j2 ) and £^(fi j2 ) 
are similar to £^(6i i2 ) and £^(6i j2 ) but with b — > r, c co / = 1 and Q — \xw) — > {\ — %w)- 



For first/second generation sfermion contributions, we find 



Ki(kn) = ^nm)J (4 5 |(T 3 - Qx w ) 



(A.20) 



where T3 is the weak isospin and Q is the electric charge assignment (taking care to flip 
the sign of Q for i?-sfermions) . For instance, 
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These contributions, arising from electroweak D-term contributions to masses, cancel out 
separately for squarks and sleptons in the limit of mass degeneracy, due to the fact that 
weak isospins and electric charges (or weak hypercharges) sum to zero in each generation. 
For this reason, we sum these contributions before taking the maximum contribution to 
the fine-tuning measure Aew- 

For chargino contributions, we find 
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For contributions from neutralinos, we find 5 
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Zi 

D{Zi) = Y[n^i( m % ~ m % )> an d TWy = M\ cos 2 Oyy+Mz sin 2 9w Our neutralino corrections 
differ in form as well as numerically from those in the literature where these were calculated 
using the neutralino mass (not mass squared) matrix [39, 40]. 
For weak bosons, we find 
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For Higgs bosons, we find 
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5 Unlike the case of other contributions where it is easy to explicitly find the eigenvalues of the Higgs-field- 
dependent squared mass matrices, this is not possible for the neutralino. To evaluate the derivatives of the 
eigenvalues of the squared neutralino mass matrix that appear in (A. 15), we use the technique introduced 
in Ref. [39] and elaborated further in Ref. [78]. 
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For SM fermions t, b and r, we find 
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